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Abstract. The main objectives of this article are two-fold. First, wc study 
the effect of the nonlinear Onsager mobility on the phase transition and on 
the well-posedness of the Cahn-Hilliard equation modeling a binary system. It 
is shown in particular that the dynamic transition is essentially independent 
of the nonlinearity of the Onsager mobility. However, the nonlinearity of the 
mobility does cause substantial technical difficulty for the well-posedness and 
for carrying out the dynamic transition analysis. For this reason, as a second 
objective, we introduce a systematic approach to deal with phase transition 
problems modeled by quasilinear partial differential equation, following the 
ideas of the dynamic transition theory developed in Ma and Wang |17l I16| . 
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1. Introduction 

The Cahn-Hilliard equation is a basic model in material science, as it character- 
izes important qualitative features of binary systems. The model has been inten- 
sively studied, especially in the case of constant mobility; see among many others 
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[11 [3 m UHl [m EU [22]. However, the dependence of the mobihty on the concentra- 
tion is very much relevant for physical applications, and a concentration dependent 
mobility appeared in the original derivation of the Cahn-Hilliard equation in [3]. 
In this case, the modeling equation is no longer a scmilincar equation, and become 
a quasilinear equation, which makes the problem much more challenging. 

The main objectives of this article arc to study the effect of the nonlinearity of 
the Onsagcr mobility on the phase transition dynamics and on the well-posedness 
of the model, and to introduce a systematic approach for studying phase transitions 
for such quasilinear systems. 

First, for a quasilinear dynamical system as the Cahn-Hilliard equation with 
the Onsager mobility, the main difficulty comes from the regularity loss through 
the nonlinear terms involving the highest order spatial derivatives. This has to be 
compensated by the regularizing properties of the linear operator. In particular, the 
so called maximal regularity property [51 [23] is essential to guarantee the existence 
of a center manifold for a quasilinear system. This can be achieved by working 
in more regular function spaces [21 [6] [TS] [20] [24]; see Section [4] for more details. 
Under this setup, wc arc able to derive the same approximation formulas for center 
manifold functions for quasilinear systems as in jl6j . With these approximations 
at our disposal, the main ideas and methods in the dynamic transition theory can 
then be applied to studying quasilinear systems. 

Second, by putting the Cahn-Hilliard equation with Onsager mobility in the 
framework just mentioned, we are able to derive the detailed transition dynamics 
as for the constant mobility case, leading to precise information on the type and 
structure of dynamic transition. In particular, we derive that as for the steady 
state bifurcation case given by Hsia , the type of transition, the critical temper- 
ature and the strength of deviation of solutions from the homogenous state are all 
independent of the choices of the nonlinearity of the Onsager mobility. 

Third, to set up the problem so that we can use the center manifold theory 
and the approximation formulas for the center manifold functions for quasilinear 
systems, we need to examine carefully the well-posedness of the model. In the con- 
stant mobility case, the equation being semilinear, the well-posedness can be dealt 
with using standard procedure for semilinear equations (see e.g. |10|). However the 
well-posedness is an issue in the non-constant mobility case and the results in this 
case are far from being satisfactory. For the two-dimensional case, the existence 
and uniqueness of a classical solution has been established recently in [T3] . But for 
the three-dimensional case, we are not aware of any such result except some partial 
results; see also [Tl[8j[23]. Hence we derive the existence and uniqueness theorems 
of global strong solution with small initial data to the equation, which is sufficient 
for the purposes of this paper. 

This article is organized as follows: The model is presented in Section 2. and 
the phase transitions for the model in a rectangular domain is given in Section 3. 
Section 4 addresses the general framework for dynamic transitions for quasilinear 
systems. Section 5 is devoted to the proofs of the phase transition results based 
on the dynamic transition theory. The existence and uniqueness of global strong 
solutions is analyzed in Section 6. 
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2. The Model 

Consider a binary system consisting of elements A and B with molar fractions 
ui and 1 — ui, respectively. The free energy of the system is given by 



aK) = ^(||Vuip + M/(rii)) dx, 



where H is an open subset in with Lipschitz boundary 9f2, a > is a constant, 
and \E'(ui), the homogeneous free energy for a mean field model of binary systems 
at a fixed temperature, is in the Hildebrand form: 

^'(ui) = RT{ui \nui + (1 - ui) ln(l - ui)) + jui{l - m). 

Here R is the molar gas constant, T is the temperature of the system measured in 
Kelvin, and 7 > is the coefficient of repulsive interaction between A and B. 

The Cahn-Hilliard equation associated with the above free energy is the follow- 
ing; see UlSllIlin]: 

^ = - V • J, 

^ ' sg 

J = - H{ui)\7fi and fi = - — = -aAui + 

OUi 

where J is the flux of type-A molecules, H{ui), a strictly positive function, is the 
Onsager mobility measuring the strength of diffusion, fi is the generalized chemical 
potential, and SQ/Sui is the variational derivative of Q. 

The above equation is supplemented with no-flux and Neumann boundary con- 
ditions: 

J ■ >^\dn =0, 
Vui • v\dn =0, 

which is equivalent to 

l^j 9n - U, — ^ — an - U, 

where v is the outward unit normal vector at the boundary dVL. As a consequence 
of the no-flux boundary condition, the mass is conserved: 

(3) ( uidx = 0. 

Now representing the deviation of concentration around a homogenous state ui 
by u = Ml — ui and approximating H{ui) and by their Taylor expansions 

about Ml, the equation governing the evolution of u can be stated as follows; see 
Hsia [n]: 



du 

at ' 

(4) - H'{ui)V [MV(aAM - biu - b2U^ + o{u^))\ 



H{ui)A [aAu - biu - 62^^ - b^u^ + o(m^)] 



iiJ"(Mi)V [u^ViaAu - biu + o(m))] 
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Here H{ui) > is the Onsager coefRcient evaluated at m = ui, and 

RT 



Ul{l - Ul) 



27, 



6. = V^ ' 



2"^ \{l^uiY u\ 

The boundary conditions in ([2]) read: 

5 ^\an = 0, ^—\dn = 0. 

ov ov 

Equation (|4]) is also supplemented with the following initial condition 

(6) m(0) = 0. 

Due to the mass conservation we assume in additional that 

(7) / itda; = 0. 

3. Effects of the Onsager Mobility on Phase Transition Dynamics 

In this section we present our theorems describing the phase transitions of Cahn- 
Hilliard equation in a rectangular box VL = nf=i(Ojii)- These theorems show the 
independence of the dynamic transition on the nonlinearity of the Onsager mobility. 

We consider the following three cases of the domain: 

(8a) L = Li> L2> £3, 

(8b) L = Li=L2> L3, 

(8c) L = Li=L2 = L3. 

The critical temperature at which the homogenous state loses its stability is given 
by: 



„ Ul (1 — Ul ) / ut/l 

(9) ^^-^Sr^ 2,-^ 



see Step 2 in Section [5] for more details. The following numbers, evaluated at Tc, 
are crucial to describe the phase transition of the problem: 

26L2 ^ 

2W2 "V ^T=T^ 



(10b) B2=[b,^^b' 

(10c) B3 = (b3 



\T=Ta 



Theorem 3.1. Assume L ^ Li > L2 > L3. Then the system (j4])-([7]) has a phase 
transition at {u,T) — (0,Tc). Moreover, the following statements are true. 
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i) If Bi > 0, then the transition is Type-I. In particular, the problem bifurcates 
on T < Tc to exactly two equilibria uf and which are attractors and can 
be expressed as 



y 3BiUi(l — Ml) L 

ii) // Bi < 0, then the transition is Type-II. In particular, the problem bi- 
furcates on T > Tc to exactly two equilibria uj and U2, which are non- 
degenerate saddle points given by: 



u 



T _^ 4R{T-Tc) TTXl 



1,2 



" 3BiUi(l — u\) L 



Theorem 3.2. Assume L = Li = L2 > L3. Then the system ([4]) -([7]) undergoes a 
phase transition at T = Tc satisfying the following properties: 

i) If B2 > 0, then the transition is Type-I and the problem bifurcates onT < Tc 
side to an attractor Yit, which is homeomorphic to the unit sphere and 
contains 8 non-degenerate singular points with 4 minimal attractors. 

ii) If B2 < 0, then the transition is Type-II and the problem bifurcates to 8 non- 
degenerate saddle points at T — Tc- There are 4 saddle points bifurcating 
out on both sides of Tc if Bi > 0, and all of the 8 bifurcated saddle points 
are on T > Tc side if Bi < . 

Theorem 3.3. Assume L = Li = L2 = L^^. There is a phase transition at 
{u,T) = (OjTc) for the system and the following assertions hold true: 

i) If B^ > 0, then the phase transition is Type-I, and the problem bifurcates on 
T < Tc side to an attractor Y^t, which is homeomorphic to the unit sphere 

. St contains 26 non-degenerate singular points, among which 

8 are minimal attractors if 63 < — at T = Tc, and 

97r^ 
22^252 

6 are minimal attractors if > 7: — at T = Tc- 

ii) If B^ < 0, then the phase transition atT ~Tc is Type-II. In particular, the 
problem bifurcates to 26 saddles at T ~ Tc. On T > Tc, there are 

8 saddle points if B2 > and B^ < , 

20 saddle points if Bi > and B2 < 0, 

26 saddle points if Bi < 0, 

and the rest are on the side when T < Tc. In all these three cases, the 
saddle points are all non- degenerate. 

Remark 3.1. When the transition is Type-II, the system undergoes a drastic 
change as T decreasingly crosses Tc. On T > Tc, the physically meaningful states 
are the homogenous state u = and some transition states away from u = which 
are metastable. The bifurcated saddles indicated in the theorems in this case are 
not physical states. 



6 



LIU, SENGUL, AND WANG 



4. Dynamic Transition Framework for Quasilinear Systems 

In this section, we present a general framework for studying phase transitions for 
quasilinear systems based on the dynamic transition theory developed recently by 
Ma and Wang [171 US] • The basic philosophy is still to search for the complete set 
of transition states as in the dynamic transition theory. For quasilinear systems, 
the key technical ingredient is the reduction of the original system to a properly 
defined center manifold for quasilinear parabolic equations [20j [24] . 

4.1. Center manifolds for quasilinear systems. Let Xi C X be two Banach 
spaces with dense and continuous inclusion. Consider 

(11) -^-L,u = G{u,X), 

U{0) = Mo, 

where u is the unknown function in C{[0, T]; X), X is a real parameter of the 
system, for each A the linear operator L\ : D{L\) = Xi ^ X \s the infinitesimal 
generator of an analytic semigroup (e^^*)t>o with domain D{L\) independent of A, 
L\ depends continuously on A, and G : Xi x M — X is a given nonlinear function, 
which contains terms of highest order derivatives in space variables and thus makes 
the problem quasilinear in nature. 

As is well known, the starting point of the existence of center manifolds is the 
variation of constants formula 

(12) u{t) = e^^'ua+ f e^^(*-")G(u(s), A) ds. 

"'0 

However, this is only a formal expression. To make sense of ([T2|) . we face two 
difficulties. First, we need the integral term to be finite and second, it should be in 
the same space as u. 

There is an easy remedy for the first one by strengthening the usual concept of 
a solution by requiring 

(13) u€Ci[0,T];X,)nC\[0,T];X). 

This requires, of course, that we choose the initial data uq in Xi. 

To overcome the second difficulty, we have to deal with the regularity loss due to 
the nonlinear term G. This has to be compensated by the regularizing properties 
of the analytic semigroup generated by the linear part. In order to achieve this, we 
have to choose our spaces carefully. As is well known (see e.g. Henry [IHI), for the 
semilinear case, this can be overcome by requiring that G : X^ x M — )• X with X^ 
being some intermediate space between Xi and X. But this does not work for the 
quasilinear case because of the terms with highest order derivatives involved in G. 
One way to fix this is to work in a pair of Banach spaces Dl^{9 + 1) and Dl^(0) 
for some 9 G (0, 1) instead of Xi and X, where Dl^{9 + 1) and D^^ (9) are defined 
as follows: 
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Definition 4.1. Let A be the infinitesimal generator of an analytic semigroup in 
X . For 6 G (0, 1), the spaces Da{()) and Da{0 + 1) are defined as: 

Da{0) = {u e X \ \\t^-''Ae''^u\\x e L°°iOA), lim \\t^-^ Ae*^u\\x ^ 0}, 

^^4^ ll"llD^(e) = \\u\\x + max^ \\t'-'> Ae'^x , 
Da{0 + l) = {ue D{A) I Au e Da{0)}, 

\\u\\DA(e+i) = \H\x + \\Au\\DA{e). 

The function spaces DA{d) and Da{0 + 1) arc Banach spaces endowed with 
corresponding norms respectively. For any 9 £ (0, 1), 

DAi0) = {X, DiA))e, 

where D{A) is the domain of A, and {X, Y)g is the real interpolation space between 
Y and X; see e.g. I5l[15l[26]. 

It is known that Da{0) does not depend explicitly on the operator A, but only 
on the domain of A and on the graph norm of A; see e.g. Corollary 2.2.3 in [15]. So 
by our assumptions on L\, Dl^{6) does not depend on A as long as A is restricted 
to some bounded interval in R. We refer readers to [T5] for some equivalent charac- 
terizations of these two spaces for arbitrary Banach space X. When X is L^'(ri) for 
some properly chosen p, these spaces are contained in the so called (little) Nikolski 
spaces hp{fl) for some s. It is this characterization and the known nice properties 
of the Nikolski spaces that help us overcome the aforementioned second difficulty. 

Now, we present the center manifold theorem for (jlip under the following as- 
sumptions: 

{Ai): The Banach space X splits into closed Lx-invariant subspaces Ei and E2 
such that (jlip takes the form 



duc 
"dT 

'df 



L^Us = P2G{Uc,Us,\), 



where u = Uc + Ug, Uc G , Ug G Xin£'2 , := ^aI^a are the restrictions 
of Lx to the corresponding invariant subspaces, and Pi : X ^ Ef are the 
canonical projections for i — 1,2. Moreover, dimi?^ < 00, all eigenvalues 
of Li have nonnegative real parts at some A = Ac, and for A sufficiently 
close to Ac the operator ■ ^1 l~l ^2 ~^ ^2 closed, densely defined and 
satisfies the resolvent estimate: 

||(L^ - z)-i£A^£A < Vz e C with Rez > 0. 

{A2): There exist neighborhoods Ui C E^ and U2 C Dj^\{9 + 1) of zero and an 
integer fc > 1 such that 

G = (PiG, P2G) G C^.,„if([/i X [/2 X M,i?i^ X D^.ie)), 

where C^^jif is the set of all functions with bounded uniformly continuous 
derivatives up to order k. Moreover, there is a neighborhood A of Ac, such 
that G(0, A) = 0, and {d/du)G{Q, A) = for ah A G A. 
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Theorem 4.1 ([2D])- Let (Ai) and {A2) be satisfied for (fTTjl . Then there exist 
neighborhoods U[ C Ui and U2 C U2 of zero, a neighborhood A' C K. of Xc, and a 
function 

with the following properties: 

i) The set 

Mx = { [u,, <^{u„ A)) e X D{L^) \ G U[], 

called the center manifold for is locally invariant, namely for each 

uo e Mx, 

ux{t,UQ)&Mx, VO<t<i„o. 
Here ux{t, uq) is the solution of (jlip with initial datum uq and t^^ is some 
positive constant depending on Uq. 

ii) $(0, A) = 0, (9/<9itc)$(0, A) = 0. 

Now we give the definitions and some crucial properties of Nifcolsfci spaces fol- 
lowing [5], from which we will see that the assumption {A2) above can be verified 
easily when we choose the spaces carefully. 

Definition 4.2 (0). Let a e (0, 1), p £ (1, 00), and n&N. Then 

/i^(R") = {ueLP(R") : \t\-''\\ui- + tej)-u{-)\\L^^O ast^O,V.? = l, ■■■ ,n}, 

where ej is the unit vector in the j*^ direction. 
For m G N and any open set f2 C K", 

h''j,{n)^{ueLP{^)\3ue /ip(M") such that u\n = u], 
- { w G I D"u G hliVl), = m}. 

Lemma 4.1 ([Sll^nj)- LetD, be an open bounded subset o/M" with smooth boundary. 

i) For s > n/p, s ^ N, the space hp{fl) is continuously embedded in C{Q) and 

thus forms an algebra. 
h) For s = m + a > n/p, m eN, a e (0, 1), and f G C""+''(R', M) with some 

fc G N, the evaluation mapping 

K(-), • • ■ , u^i-)) G ih;in)y ^ • ■ • , e /i^(fi) 

zs times continuously differentiable. 

We note that the first part of our assumption {A2) is a direct consequence of 
Lemma 14.11 under the condition that G is smooth enough. 

4.2. Approximation of the center manifold function. In this subsection, we 
consider an approximation of the center manifold function A) for (|lip obtained 
by Theorem 14. II following the same line as in [TB]. 

We assume that the nonlinear term G{u, A) in pT|) has the Taylor expansion 
about u = as follows 

+ o(||u||^^ (9+1))' foi" some 2 < k < r, 

m—k 

where u G Dl^{0 + 1), G,„ : L»l^(6I + 1) x ••• x Dl^{6 + \) Dl^{0) is an m- 

" V ' 

m times 

multiple linear operator, and Gm{u, A) = Gm{u, • • • ,m, A). 
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Let {/3i{\) e C I 2 e N} be all eigenvalues of L\ counting multiplicities and 
{ei(X) I z G N} be the corresponding eigenvectors. Assume that the following prin- 
ciple of exchange of stabilities (PES) condition holds: 



(16) 



< if A < Ac 
Re^»(A)<(=0 if A = Ac Vl<i<m, 

> if A > Ac 
Re;3j(Ac)<0 Vj>m+1, 

for some Ac G M. 

Wc also assume that the span of {ei(A) | i G N} is dense in DL^{d + 1); namely 

(17) DlAO + 1) - span{e,{X) \ i G 
Now, let 

= span{ei(A), • • • , em(A)}, 
E2 = the complement of E^ in X. 
Then Lx is invariant on E^ and E2, i.e., L\ can be decomposed as 

L\ ~ Li (B L2, 

(18) : E^ ^ E^ 

L2 ■■ Xif] E2 E2, 

where is the Jordan matrix of L\ associated with /3j(A) {1 < i < m), and 
has eigenvalues /3j(A) {j > m + 1). 

Now, we present the following theorem which gives a first order approximation 
formula of the center manifold function of (|lip for A close to Ac . The approximation 
formula is essential to understand the dynamic behavior of the trivial solution u = 
of PT|) for A near Ac. 

Theorem 4.2. Assume all the above conditions given in this subsection hold. For 
the nonlinear term G(u,X), assume in addition that {A2) in Subsection \4.1\ holds. 
Then for A sufficiently close to Ac we have the following approximation for the 
center manifold function $(uc. A).' 

rO 



(19) ^u,^ A) = / e-^^^P2Gfe(e"^i«c, A) dr + o(||«c|l?,, je+i)), 

J —00 

where and are the linear operators as given in (jlSp . Gk{u, A) is the lowest 
order k-multiple linear operator as in (|15|) . and = X^I^li 2/*^* ^ ^Lxi.^ + 1) 
is sufficiently .small. In particular, for some special cases we have the following 
assertions: 

i) if Li is diagonal near A — Ac, then (jl9p can be approximated as 

(20) -L^^u,,X)^P2Gk{uc,X) + oik). 
Henceforth, o{k) stands for 

(21) o{k) := o(|kc||^,^(,+i)) + 0{\ReP{X)\\\uctD,^^e+i)), 
with /3(A) being the eigenvalue of Lx with largest real part. 
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ii) Let m = 2 and /Ji(A) = = a(A) + ip{\) with p{\c) ^ 0. // Gk{u, X) 
is bilinear, i.e. k = 2, then the center manifold function $(mc, A) can he 
expressed as 

[(-L^f +V(A)](-i^)$K, A) 

= [{-L^f + V(A)] P2G2K, A) - 2p\\)P2G2{u,, A) 
(22) +2p\\)P2G2{yie2~V2ei,\) 

+ p(A)(-L2)[P2G2(yiei + ^262, 2/261 - yie2, A) 
+ P2G2{y2ei - 2/162, 2/iei + 2/262, A)] + o(2). 

iii) Let /3(A) = /3i(A) = /3m(A) /lawe algebraic multiplicity m > 2 and 
geometric multiplicity r ~ 1 near A = Ac, i.e., has the Jordan form: 

( f3{\) 5 ••• \ 
/3(A) ••• 



(23) 



Li 



Let 





V 



/3(A) 5 

/3(A) y 



for some 5^0. 



j=l r=0 

where y = (i/i, ■ • • , 2/m) G 1^'", S is as in (|23p . and t > 0. T/ien t/iere exist 
functions i^o(2/), ' ' ' , Fk(m-i){y) such that the k-linear term Gk{z, A) can 
be expressed as 

Gfe(z, A) = Fo(2/) + iFi(2/) + • • • + J^fe(™-i) (2/), 

and f/ie center manifold function $ /las i/ie following form 

fc(m-l) 

- (i^)^+i<i>, = j!P2Fj(2/), for 0<j< fc(m - 1). 

The above Theorem is a direct generahzation of the Hilbertian version in |16| 
and the proof is the same as the Hilbertian version with obvious modification and 
is thus omitted here. 

5. Proof of Main Theorems on Phase Transitions 



(24) 



In this section, we provide a unified proof for Theorem I3.1H3.3I on the phase 
transitions of the problem (l4])-(l7]). The main ingredient of our proof is the center 
manifold reduction, following the line of Ma and Wang [TH] . But since our equation 
is quasilinear, it seems very hard, if not impossible, to do the reduction in Hilbert 
space setting as was done for semilinear case in [TH]; see also the discussion in 
Sectional Instead, we will work with a pair of Banach spaces (7?^^ (0 + 1), (6*)) 
for some 9 G (0,1) as defined in Definition 14.11 where the existence of a center 
manifold is known and is recalled in Theorem 14. II 

In order to study the phase transition of the problem we need that the equation 
admits a global solution u £ Cb{[0, 00); Dl^{9 + 1)) n C^([0, 00); Dlt{0)) at least 
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for small initial data in DLj,{d + 1). This is done in Section IH where the existence 
of global solutions with small initial data in is also shown. 

Assuming for the moment the wcU-posedness of the problem (H))-© with small 
initial data in (0+1). we prove the main theorems in five steps. In the first step, 
we establish the necessary functional set-up. In Step 2, we analyze the linearized 
problem to identify the critical parameter at which the homogeneous state w = of 
the system loses its stability. Step 3 is devoted to deriving an approximation of the 
center manifold function by the approximation formula given in Section 14.21 We 
derive the reduced equations to center manifolds in Step 4. In the last step, the 
reduced equation to the corresponding center manifold is analyzed. 

Step 1: Functional setting. For the functional setting of the problem, we will 
choose p > 3 and 6 > such that 1 > 49 > 3/p and set 



(25) 



D{Lt) ^{ue W^'^m \ ^^^ = of udx = 0}, 

X = {ue LP{n) I / udx = 0}. 
Jn 



With this choice of p and 9, the interpolation space Dlt{0) in (ESI) becomes an 
algebra (see Lemma 231, which is essential to guarantee the existence of a center 
manifold. We note that the algebra property is also needed for the well-posedness; 
see the proof of Theorem 16.21 

We define the operators Lt = —A + Bt ■ D{Lt) X hy 

Au = aH(ui)A'^u, 
^ ' Btu = biH{ui)Au, 

and G by 

G{u, T) ^H{ui)A{b2U^ + b3U^ + o{u^)) 
(-27) - H'{ui)V [uV{aAu - biu - &2U^ + o(u^))] 

- ii7"(ui)V [u^V{aAu - bm + o{u))\ . 

The problem jlJ-Q can now be recast in the following abstract form: 
du 

(28) —=Ltu + G{u,T), u{0) = ^. 

at 

Letting s = 49, it is known (see [5|) that the interpolation spaces Dl^{9) and 
Dlt{^ + 1) defined in Definition 14.11 are given by 



DLA&)^{X,D{LT))e^{u^hl\ [ udx=^0}, 

Jn 

DlAo + i) = {ue h;+\n) \ — = ^ = o, udx^o}, 



(29) 



where is the Nikolski space defined in Definition [ 

From Lemma 14.11 we know that for s = m + a > n/p, < cr < 1, / G 
C"'+'^(R', K), the evaluation mapping 

(^i(.), • • ■ , ui{-)) G {h;{n)y ^ /(ui(.), • ■ • , ui{-)) G h;{n) 

is k-times continuously differentiable. This immediately implies that 
(30) G(-,r) : Dlt{9 + 1) ^ Dlt{9) is smooth for all T > 0. 
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We can also check easily that 

(31) G(0, T) = and — G(0, T) = for aU T > 0. 

ou 

Step 2: The principle of exchange of stabilities (PES). In this step, we explore 
the eigenvalue problem associated with the linearized counterpart of psp to identify 
the critical parameter T = Tc at which the homogeneous state u = of the system 
loses its stability. First, we consider the eigenvalue problem 

- AcK = PkEk in 

dcK. 

(32) ^1^"-^' 

ck da; = 0. 



The eigenvectors ck and the eigenvalues pK are given by 

(33) Cif = J|cos^— ^, PK = ^-jjr^ 

z=l ' i=l ^ 

where 

KelC:^ {(fci, fc2, fca) : h >0,kf + k^ + kl ^ 0}. 
Now, we turn to the eigenvalue problem associated with the linearization of 
around it = 0: 

(34) LtCk - pK{T)eK. 

It is easy to see that the eigenvectors of (p4| are the same as the eigenvectors of 
([32]) . and the eigenvalues are given by 

= H(ui)pK 27 - — — - apK 

\ Ul(l-Ul) 

Let Tc be given by One can readily see that (3k (T) < for all K G K, when 
T < Tc. Now, we define V, a subset of JC, which contains all K G K. satisfying 
/3if (Tc) = 0; namely 

'{(1,0,0)} ifLi>L2>T3, 
{(1,0,0), (0,1,0)} ifLi = L2>T3, 

[{(1,0,0), (0,1,0), (0,0,1)} if Li = L2 = ^3- 
By (|33|) . ([35]) and our choices of Tc and V, we see that PES is valid: 

r < if T > Tc, 
(3k{T)}^0 if T = Tc, yKeV, 
[>0 if T<Tc, 

/3x(Tc)<0, yKeic\v. 

The PES above shows that Tc is the critical parameter value at which the homoge- 
neous state loses its linear stability. From the general dynamic transition in [16] , we 
know then that the system will always undergo a dynamic transition at this critical 
threshold. The type of transitions is however dictated by the nonlinear interactions, 
which we shall explore in the next few steps. 



(36) V : 



(37) 
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Step 3: Approximation of the center manifold function. Let E-f = spanjcif | 
K G V} and E'J be the complement of Ef in X, where X is defined in ((25|). Let 
and L2 be the restrictions of Lt to Ef and respectively. It is clear that 
assumption (^1) below Definition 14.11 is satisfied for ((28)) with T playing the role 
of A. Thanks to ^ and (|3T]), {A2) is also satisfied. Thus, by Theorem O the 
system ([28]) admits a center manifold in a neighborhood of u ~ in Dl^{6 + 1). 
In the following, we will use Theorem 14.21 to derive an approximation of the center 
manifold function $(uc,r). 

Let 

(38) Uc=Y. V'^J^ 

.lev 

where = Vivi^yi^ for J = (.71, .72, .73)- Let 

(39) S = {J + L\J,LeV}. 

For example liV = {(1, 0, 0), (0, 1, 0)} then S = {(2,0,0), (1,1,0), (0,2,0)}. 

Note that the Jordan matrix if is diagonal for all the three types of domain 
n as given in (|8a| -([8c|). then we have the following approximation of the center 
manifold function $ (see Section |4?2] formula (f20| ): 

(40) -L^$(u„r) = P2G2(uc,r) + o(2), 
where G2 consists of the quadratic terms of G given in ((27)) . i.e., 

(41) G2(uc, T) = H{ui)b2Aul - iJ'(?Ii)V(wcV(aAuc - hua)), 

and the notation o{n) is as in (PT|) with T playing the role of A. Henceforth, all the 
equalities involving T hold for T sufficiently close to Tc- 

Let (•, •) denote the L2 inner product. Note that we have the following orthogo- 
nality relations 

V 0, if J = K, 
= 0, if J ^ K. 

Since <&(mc, T) £ i?J and {bk \ K G IC\P} spans iJj, by the orthogonality relations 
above, we can write $ in the following form: 

Now, for each ck with K € JCXV, we take the L2 inner product of PO]) with e/^ 
and integrate by parts on the left hand side to obtain 

(-i^$(u„T),eK) = -mu,,T),L^eK) 

^-l3Kmu„T),eK) = {P2G2{uc,T),eK)+o{2) 

- {G2iuc,T),eK)+oi2). 

The last equality above holds due to (H^ . Thus, 

\ {G2{uc,T),eK) 

PK 

Plugging this back to we obtain 

(44) 4>(..,T).-Ei%(^^l^e. + o(2). 



(42) (ej,eK) 
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Also note that for all K e V, Pk(T) ^ as T ^ Tc. Then by §^ and ([35]), we 
have 

(45) apK + &i = 0{13k{T)) as T ^ for all K e V. 
We now compute the term {G2{uc,T),eK)- By ([38]) . we have 

(46) {Aul,eK) = {ul^Aex) = -PK ^ / y^^^ cjClCk Ax, 
and 

(V • (ucV(aAwc - huc)), ck) 

= -iY. y'^J^iYl /(aAei -6ieL)),Veif) 
Jev Lev 



(47) = y'ejViJ2 + &i)eL), Ve^) 

Je-p Lev 

= ^ {apL + bi)y-''^^ / ejVeLVe^ dx. 

By our definitions of V and <S in (p6|) and (|39p , respectively, one can easily see that 
for any given J, L E V and A' £ /C \ "P, we have: 



(48) 



ejeLe^da; 



n 



/ ejVeiVe/f da; 



I ifK = J + L, 

otherwise, 

f^^O if/v=J + L, 
= otherwise. 



Here V — L1L2L3 is the volume of ft. 
Now, by dH]), and (l46l)-(|48l), we have 

(49) (G2(wc,T),eK) = 0, V/v e /C\ (T'U^). 
By (gS]) and (gT]), we also have 

(50) (V • (iieV(aAu, - hu,)), ck) = o(2), V if e 5. 

Hence by pS)) . (|^^ and ([SO)) , the center manifold has the following approxi- 

mation: 

Hy) = ^KeK + o{2), 
Kes 

(^1) (g(lIi)fe2A^g,eK) H{u,)b2PK ^ /• H ^ <^ 

= in ^ — in \ 1^ y ejeLCRdx, K eS. 

-pK{eK,eK) PK{eK,eK) j'j^^ Jn 

Using dSS]), (gH) and gll), we have 

1>2,7 = -^ + 0(/3i(r)|2/2./|), Jep, 

$j+L = + 0(/3i(r)|y^+^|), J^LandJ, Lev. 

apj 

Step 4: Derivation of the reduced system. Now let 
(53) M= ^/e,7 + <f>(2/,T). 
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The dynamics of the system ([28| close to Tc is determined by the dynamics on 
the corresponding center manifold. To this end, we replace u in ([28]) by the right 
hand side of ([551) . take the L2 inner product of ([^5]) with ej, and make use of the 
orthogonality relations (|42p to obtain the following reduced system: 



(54) 



at {ej,ej) 



J eV. 



The second term on the right hand side of (|54p can be simplified further using the 
approximation formula of the center manifold function (|52p as we now show. For 
J P, making use of the orthogonality relations (|42p . the following can be obtained 
by direct computation: 



(55) 



(56) 



(57) 



(58) 



(Au^,ej) = -/9j(u^,ej) = -2pj y'^^K / eLCRej dx + o{3) 



Lev,Kes 
fj+L -I- 0(3). 



L6P 



(A^i^eJ) = -p., 



-PJ 



Wp 



K,L,M£V ''^^ 

{y'''fejdx + 3 2/^+'^ /44dx)+o(3) 



LeV,L^J 
J+2L 



y 



0(3). 



L£V,L^J 



(V • {u\7{aAu - hiu),ej) 

= Y y^^K{apK + hi) I eiVexVej da; + 0(3) 

V y^^j+L{ap.j+L + &i) / eLVe,/+L Ve,7 dx + o(3) 



[2/(ap2j + 6i)*2j + Y y^'^J+dapj+L+hi] 

LeV.L^J 



0(3) 



(V-(?/Vu2),e,;) 



E 



K+L+M 



K,L,MeV 



/ e^eiVeM ■ Vejdx+ (u^, Aej) +0(3) 



,J+2L 



|Vej|d.T + (u^Aej) +0(3) 



(y3^ + 2 ^ /+2^) + (u3,Ae,,)+o(3) 



(2/3.^ + 2 ^ /+2^) 



0(3). 
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The last equality above follows from the result of ([SG]). 

(V • {u'^V{aAu - biu)), ej) 

= -( E y''^''^KeLV{ J2 /'(aAeM - ^iCa/)), Vej) + o(3) 
K,Lev Mev 

(59) 



J2 y'^+'^+^'iapM + bi) j CKeLVeMVej dx + o(3) 



K,L,M£V 

o(3) (by (113 



Using (I5T])-(I52]) and (l55])-(|59l) in <^ and p.j = for all J e P with L as in 
(|8ap ~ ((5c) ). we get the following reduced system: 

(60) ^ = PjiT)y'^^^^y'ia^y''^ + a,J2y'^) + oi3), JeV, 

Lev 

where ai and (T2 are 

ID^i f^i 5-; and (72 = 363 5-- 

Step 5: Analysis of the reduced system. The reduced equation (|60p is essentially 
the same as in the case of constant mobility except for a factor of H{ui) appearing 
in the cubic terms; see Ma and Wang [l8l. For the sake of completeness, we present 
here the main ingredients of the analysis. 

First, it is known that the transition type of (|60| at the critical point Tc given 
by ^ is completely determined by the following equations: 

(62) '^--^^yH4y'' + 4Y.y'') vJeP, 



Lev 
L^J 



dt 2L 
where 

(63) crj = (Ti|t=t<,, and (T2 ~ a-2\T=Ta- 

Recall Bi, B2 and given in (|10ap - (|10cp . It is easy to see that 
(T? > ^ Bi > 0, cr? < ^ Si < 0, 

(64) Cr? + (72 > ^ B2 > 0, (7? + (7^ < ^ B2 < 0, 

(7? + 2(7^ > ^ B3 > 0, crj + 2(7^ < ^ B3 < 0. 

These relations will be used frequently in the following. 

Second, for the case where L = Li > L2 > L3. the critical index set P = 
{(1,0,0)}, the equation reads: 

(65) = /3(i,o,o)yi 5^2— ^1^1 + "(3)' 

and ([5^ takes the following form: 

c^yi H{ui)n^ 3 

Thus, the system has a pitchfork bifurcation at Tc, and the type of transition 
depends on the sign of cr?. If > 0, namely Bi > 0, then the bifurcation happens 



CAHN-HILLIARD EQUATION WITH ONSAGER MOBILITY 



17 



on the side when T < Tc, the bifurcated two steady states are local attractors, and 
the transition is Type-I. If cr^* < 0, namely Bi < 0, the bifurcation happens on the 
side when T > the bifurcated two steady states are both saddle points, and the 
transition is Type-II. It is clear now that the assertions in Theorem 13.11 hold true. 



Third, for the case where L ^ Li = L2 > L^, the equations in (|60|) read: 
-TT = /3(i,o,o) (7^)2/1 ^1^— yi(f^il/i + '^22/2) + o(3). 



^ lit- \ 2 

= P(QA,Q){T)y2 ^^5— ?/2(cri2/2 + cr22/i) + 0(3), 



and the equations in ([62]) read: 



dyi _ H{ui)tt^ 2 



(68) 



dt 2L2 



yiKyf + ^2y2), 



-y2[(JiV2+<^2yi)- 



dt 2L2 

To analyze (j68p . we first find the straight line orbits, which arc orbits of the form 

2/2 = TO12/1 or 2/1 = "122/2- 
We assume that the line 

2/2 = miyi 

is a straight line orbit of ((68|) with some mi G M. Then 

^fto^ '^y^ '^""^1 + ^^2 

(69) — = mi = jTii-fT— — — 2. 



dy 



Thus mi = 0, ±1 provided (Ti 7^ ctj- Similarly, in order that yi — m2j/2 be a 
straight line orbit of 7712 can only take the values 0, ±1 provided 7^ cr^. 

There are four straight lines in total determined by ?/2 = "^12/1 and yi = m22/2 
with mi, m2 = 0, ±1, and each of them contains two orbits. Hence, the system 
([68]) has exactly eight straight line orbits provided that CTi 7^ cr^. 

Since (|67p is a gradient- type equation, the energy decreases along the orbits. 
Therefore there are no elliptic regions at ?/ = 0. Hence, when cti + > and 
(7^ 7^ f 2 all the straight line orbits tend to ?/ = which implies that the regions are 
parabolic and stable, therefore y = is asymptotically stable for ((67|) . Accordingly, 
by the attractor bifurcation theorem. Theorem 6.1 in [17], the transition of ([67]) at 
Tc is Type-I. 

When (Ti = 0-2, one can check directly that (Ti = > 0. In this case, it is clear 
that ?/ = is an asymptotically stable singular point of ([55]) . Hence, the transition 
of ([67]) at Tc is Type-I. 

When CTi + (T2 < and (7^ > 0, namely i?i < and B2 > 0, the four straight 
line orbits on y2 = ±?/i extend outward from y = 0, and the other four on ?/i = 
or 2/2 = go toward y ~ which implies that all regions at y ~ are hyperbolic. 
Hence, by Theorem A.3 in [TB], the transition of ([S7]) at Tc is Type-II. 

When (7° ^ 0, then CTj < too. In this case, no orbits of ([S5]) go toward y ~ 
which implies by Theorem ^.3 in [18] that the transition is Type-II. 

Thus by ([64]) and the above analysis, we proved that the transition of ([70]) from 
{u, T) = (0, Tc) is Type-I if B2 > 0, and Type-II if B2 < 0. This proves the 
assertions about the types of transitions stated in Theorem 13. 21 
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Fourth, for the case where L = Li = L2 = L-s, the equations in ([60]) read 

2/i[^?y? + '^2(yi + + 0(3), 
y3[<^bl + <^2{yl + yi)] + o{3), 

and (p^ are written as 

yiW°iy! + <^"2iy'2 + yl)], 
y2[<^W2 + 4iyl + yl)l 
y3Wiyl + <^2{y! + yl)]- 

It is clear that the straight hnes 

yi = 0, yj = 

(72) yf = yl yk=0 

222 

2/1=2/2= 2/3' 

consist of orbits of ([7T|) . There are 13 straight hnes in total contained in ([7^ . each 
of which consists of two orbits. Thus, ([7T|) has at least 26 straight line orbits. In 
fact, as shown in [18], the number of straight line orbits of ([71]) is exactly 26 when 
a? ^ ct". 

As before, when ai = (J2, we have that (Ji~a2> 0. In this case, it is clear that 
?; = is an asymptotically stable singular point of (j7ip . Hence, the transition of 
([7(7| at Tc is Type-I. 

When a I + > and o^i (J2, all straight line orbits of ([71]) go toward j/ = 0, 
which implies that the regions at y = 0, arc stable, and y ~ is asymptotically 
stable. Thereby the transition of ([70]) is Typc-I. 

When (Ji + 2(72 < ^-^^^ ""i + ^2 > we can check that ^ ctj ^^nd hence all 
straight line orbits of ([TT]) are given by ([7^ . Moreover, all the straight line orbits 
determined by 2/1=2/2= 2/1 extend outward the origin, and all the rest straight 
line orbits go toward the origin. Hence, for any initial data in a small neighborhood 
of 0, the orbit of ([71]) goes away from as long as the initial data does not belong 
to any of the coordinate planes, which implies that the transition is Type-II. 

Similarly, when trj + CTj ^ one can also check that given a small neighborhood 
of 0, there is a dense subset of the neighborhood, such that for any initial data 
in the dense subset, the orbit of ([TT]) goes away from 0. Hence, the transition is 
Type-II. 

Thus by ([64]) and the above analysis we proved that the transition of ([70]) from 
(u, T) = (0, Tc) is Type-I if B3 > 0, and Type-II if B3 < 0. This proves the 
assertions about the types of transitions stated in Theorem 13.31 

Fifth, we show the nondegeneracy of bifurcated steady states. Since the bifur- 
cated equilibrium points of (j28p are in one-to-one correspondence to the bifurcated 
equilibrium points of (|60p. it is sufficient to consider the leading order steady state 



ayi n irr-. 

= P(i,o,Q) (7^)2/1 - 

(70) ^ = /3(,^^,)(r)y2- 

-77 = P(0,0,l)(^)2/3 - 



(71) 



dy\ 
dt 
dy2 
dt 
dys 
dt 



for i ^ j, 1 < i, j < 3, 

for i / j, i ^k, j ^k, 1 < i,j, fc < 3, 
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equations of the reduced system (|60| 

(73) /3j(r)y'^-/(aiy2'^ + a2^?/2^) = for J e T', 

Lev 

where aj = H{ui)TT'^iJi/{2L'^),a2 = H{ui)n'^a2/{2L'^). 

Let m ~ \V\. In [T5], it is shown that ([75)1 has 3'" — 1 bifurcated solutions, and 
all bifurcated solutions of ((73|) are regular. 

For Typc-I transition case, since all bifurcated singular points of (|28|) are non- 
degenerate and when St is restricted to j/ij/j-plane (1 < i,j < to) the singular 
points are connected by their stable and unstable manifolds, all singular points in 
St are connected by their stable and unstable manifolds. Therefore, St must be 
homeomorphic to a sphere S"™"^. 

Finally, in addition, as in |18| . the number of minimal attractors is obtained 
by studying the Jacobian matrix of (f73| . The proofs of Theorems I3.1H3.3I are now 
complete. 

6. Existence and Uniqueness of Global Strong Solutions 

In this section, we will give two results concerning the existence and uniqueness 
of solutions with small initial data, one in Hilbert space setting and the other in 
the interpolation space setting. 

6.1. Existence in Hilbert spaces. We start with the following problem: 

— = V • [i7(u)V(-aAu + hiu + b^u'^ + 63^^)] 

du dAu 
-^9^ = 0, -- — ao = 0, 

(74) ov av 

u{0) = uo. 



u 



Ax = 0. 



Here a, 6i,&2, and 63 are constants with a > and 63 > 0, and £7 is a bounded 

domain in 'S? with sufficient smooth boundary. 

We make the following assumption on the Onsager mobility H{s): 
{H): min H{s) > i?i > 0, and H{s) and H'{s) satisfy the following growth 
condition: 

\H{s)\ < C{\s\P+^ + 1), \H'{s)\ < C{\s\P + 1) V s G R, 
where 1 < p < 3. 

It is clear that the free energy functional associated with ([7^ takes the following 
form (see Section [2]): 

(75) G{u)^ f ^\Vu\'' + hiu^ + h2U^ + ]b3U^dx, 



_ 2' ' 2 3 4 

and the generalized chemical potential /i in this case is given by 

(76) := — = -aAu + biu + b2u'^ +b3U^. 

ou 

We use the following notations. | • | denotes either the norm on L^{il) or the 
Euclidean norm on R", which should be clear from the context, \ ■ \x denotes the 
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norm on the generic Banach space X, (•, •) is the L^(ri) inner product, -ff™ is the 
usual Sobolev space, and we also denote: 

(77) := {w e i/^(l^) I ^lao = and [wdx = 0}, 

Jn 

(78) M^i :={it'efl"*(f^)|^|ao-^|ao = Oand fwdx^O}, 

(79) A{t):^\u{t)\]j^+l. 

Hereafter, C denotes a generic constant which depends only on the bound Bi, 
the coefhcients foi, 62, ^3, and the domain $7, C(uo) denotes a generic constant 
depending on the initial data uq. 

We have the following existence and uniqueness theorem of a strong solution to 
the problem ([74| . which will be proved in Section [631 



Theorem 6.1. There exists a constant eq > 0, such that for any initial datum 
uq with |ito|_f/2 < eo, there exists a unique strong solution u to (|74p such that 



u e 2.2(0, T; M^i) nC([0, T]; W) with ^ e L'^{0,T; L^{n)) VT > 0. 

6.2. Existence in interpolation spaces. Now recall the Cahn-Hilliard equation 
with Onsagcr mobility: 

(80) 'i-Lru^G(u.n 

u(0) = Uq, 

where Lt is as in p6|) and G is as in p7|) . The main result for (|80)) is as follows: 



Theorem 6.2. Let Dl^{6) and Dl^{6 + 1) he as in with some p > i and 

6 > such that 1 > 49 > 3 /p. Then 3 e > and r > such that ^ T > Tc — e, 
V Uq G -8(0, r) C Dl^{9 + 1), the equation (|80p has a unique strong solution u £ 
C6([0, 00); DL^i9 + 1)) n Ci([0, 00); DL^{e)) with u{0) = uo- 



The proof of this theorem will be given in Section [6.41 

6.3. Proof of Theorem 16.11 The proof is carried out by first proving a local 
existence result. For this purpose, we need the following lemmas. 

Lemma 6.1. |Am| is a norm on W which is equivalent to the H^-norm. Similarly, 
A^uj is a norm on Wi which is equivalent to the H'^-norm. Moreover, for any 
u £ Wi, there exists a constant C depending only on the domain such that 

(81) \u\h^ < C|VAm|. 

Proof. The above results follow from the regularity theory for elliptic boundary- 
value problems. For the first claim, we use the regularity theory of the Neumann 
problem 

du 

Au = hmn, -r-|so = 0, 
on 

which implies that 

\u-iu)n\H- <C\h\^C\Au\, 
where (u)r2 = '\n\ In^'^^' ^i^ce each u €: W satisfies J^^udx — 0, the first claim 
follows. 
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The second claim follows from the regularity theory of the Neumann biharmonic 
problem 

A u = hmil, = -g;^\dn = 0, 

which implies that 

\u-{u)n\H^ <C\h\ ^C\A'^u\. 

For more details, we refer the interested readers to [25], Chapter III Lemma 4.2. 

For the third claim we use a special case of Corollary 27 in [7], which states that 
if u e H^{n) and f^lao = 0, then there exists a constant C depending only on Q 
such that 

\u\hs < C\Au\hi. 

Thus, for any u £ Wi, 

\u\h^ < C\Au\hi < C(|VAu| + |Au|) < C|VAu|. 

The last inequality follows by applying the Poincare's inequality to Am and making 
use of the fact that Au has mean zero due to Gauss divergence theorem and §^ |ao = 
0. 

□ 

Lemma 6.2. Let u{t) be a solution to (j74p with initial data uq G W. Then we 
have the following estimates: 

(82) giuit))<C\uo\ji2il + \uo\l2), Vt>0, 

(83) \uit)\m <Cil + \uo\l2), yt>0, 



(84) 



^ \u{T)\j,sdT<C\uofH.{l + \uo\h) 



eC{l + \uo\H2y°, Vt > andO < e < 1. 



Proof. Taking the time derivative of the free energy functional given in (|75p and 
using assumption (H), we have 

^ =(S' "*) = (M' V • iHiu)V^l)) = - / H{u)\Vfx\^ dx < 0, 



where fi is as in (|76p . Thus, 

GHt)) <g(u(o)) 



^^^^ <^\uo\l. + \uo\l^ f {b,ul + C)dx 

<C|uo|H2(l + |uo|i2), yt>o, 



which justifies ([82]) . 

By dTS]) and ([SS]), we have 

/ (f IV^I' + l^iu^ + 1^2^' + \b3u' ] dx = g{u{t)) <C\uo{i,2il + l^oliO, 
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which implies 

a\Vuf + h3\u\i,< I (|6i|u2 + ^|62||«|')dx + C|«o|H2(l + |«o|i2) 

<jjy{^^ + \hu^) + (c| + \bM')) dx + C|«o|l,.(l + l^oliO 

Wc thus obtain 

|V7/|2 < C\m\]j2{l + + C<C{1 + \uo\]j.f, 



and ([85)) follows by the Poincarc's inequality. 

Recall that /i = — aAu + 6iu + b2U^ + 63?/'^. We have by triangle inequality 

a|VAup < 2|V^P + 2 / |V(6iu + &2U^ + foait^)^ dx 
<2|V/^|2 + C(|Vu|2 + \u\l, \Vu\l, + Iwlts \Vu\l,) 
<2|V^P + C|u|^i + C\u\]j, |VAu|3 + C|u|S |VAit|3 
<2|VAip + C|it|^i + C\u\% |VAu|3 + C|u||i |VAu|3 |VAw|3 
<2|V/i|2 + C\u\]j^ + C\u\% + ||VAu|2 + C\u\^^i + ||VAu|2, 

where in the second last inequality we used the interpolation inequality ImI"* 9 < 

He. 

JJ3 1 

C|u|^i |u|^3 and the fact that \u\[ja is equivalent to |VAw| as shown in Lemma |6. II 
Then 

(86) I VAu|2 <C\V^l\^ + C{\u\^^i + 1). 

Note also 

pt + e n nt + l 



Bi / \Vi^ifdxdT< / / H{u)\\7fj,fdxdT 

Jt Jn Jt Jn 

))-g{u{t + 

<C\uo\l.a + \uo\l.). 



^^^^ =^g{u{t))-g{u{t + i)) 



By dMl), (EH) and dHT]), we have 

ft+e 

|VA?/|2 dr <C|mo|h2(1 + |uo|i2) + eC(l + |uo|ff2)i°. 
Now ([H^ follows from ([55)) and the fact that |VAm| is an equivalent norm to 

□ 

With the above two lemmas at our disposal, we are ready to prove the following 
local well-posedness result. 
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Proposition 6.3. For any initial datum uq G W , there exist Tq > and a unique 
local solution u(t) to the problem (|74p such that: 

(89) 

u e L^O, To; Wi) n C([0, Tq]; W) with ^ e {0 , To; L^n)) , 

at 

^{u, v) ^ I {H{u)iiAv + H'{u)^iVu ■ Vw) Ax, y v e W and a.e. < t < To, 
m(0) = uo- 

Proof. The proof consists of several steps. 
Step 1 . Given any m S N, let 

(90) t¥™ = span{efc I 1 < fc < m} C W^n ^ C\[0 , T^], Wm) , 

where e^'s are eigenvectors of —A with Neumann boundary condition on dVl and 
Cfc da; = 0, and > is a constant to be chosen as follows. 
According to standard existence theory for ordinary differential equations, for 

each m, there exist Tm > and an approximate solution u„i to (|89p in the following 

sense: 

m 

Um = ^ X]{t)ej e W,n, Xj{t) G M, 
3=0 

d f 

(91) -TjiUrn, w) ^ I {H{Um)^ynAw + H' (Um) f^m'^^Um ■ Vw) dx, Vw £ Wm, 

m 

'Um(O) = ^{uo, ej)ej, 

where flm = -aAUm + biU„i + &2Wm + &3"m- 

In order to show that there exists a solution to the original system, we need 
to establish some uniform estimates on the approximate solutions, which is the 
direction that we turn now. 

In (|9ip . using A^w,„ as the test function, integration by parts twice and applying 
(H), we obtain 

i^|Au,„|2 + aBi|A2lt„,p < {H{Um,)A{biUra + b^uf^ + bsul,), A^w™) 

■■=h+l2- 

We have the following estimates for Ii and l2- 

h ^{H{Um)A{biUyn + &2Wm + hu%^) , A^W^) 

(93) /■ 

<C {\u^\P+^ + l)\A{b^u„, + b2ul^+b^ul^)A^uJdx. 
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Note that 

= / IhilUmf'^^ + l){3ul^Au,n + 6u,n\\7u,nf)A'^U,J dx 

Jn 

<C{\Um\L~^ + l)|Um|i=o|Aw„J|A^W„i 

+ C{\umfj^ + l)|w,„|L=o|Vwm|^4|A^Um| 
<C{\u,^f+^ + 1)\u„,\UA^u,„] 

Similarly, we have 

/ \bi{\u^\P+' + l)Au,nA^u^\ dx < ^\A^u^\^ + + lf\u. 

Jn 



2 



and 

/ HlUml^^' + l)Aul,A^Um\ dx < ^\A\„-,\' + C i\UraC,' + l)^\u^\%.. 

Plugging the above three inequahties into (1931) . wc have 

<^\A\^\' + ci\u,^\l.. + ir+\ 

For I2 , we have 

h ={H'{Ur,i)VUr,i ■ V/im, A'^Urn) 

= {-aH' {Ujn)VUjn ■ VAUto, A^Ura) 

+ (i/'(lt„)Vu„j • V{biU,n + b2u'^a + A^ltm). 

By our assumption on H, the first part of I2 can be estimated as 
(-aiJ'(u,„)Vu„,-VAu,„, A^Uryi) 

<C{\u^\l^ + 1)|Vu™|loo|VAw,„||A2u„,| 

<(by Agmon's inequality, see e.g. [IS] page 52) 

1 1 

Umffj2 + 1)|Vm,„||ji |VUrn||f2 |Um|//3 |A Um\ 

- - 9 

<C{\Ura\ ^2 + 1) I jj2 I jjs I A Uyn\ 

<C(|u™|Pj2 + l)\u^\%^\A^u„,\i 
The second part of I2 can be estimated in the same fashion, and we have 
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By combining the estimates for the two parts of I2 , we obtain 
'4 



h <^^|A^W,„P + C(|u„i|^2 + l)^(|Wm||f2 + \u,-n\%2 + |Um||f2) 



(95) +C(|u™|^. 



By ([921), ([94]), ([951) and Lemma 16711 we have 

(96) ur 
Set 

2/ = 1 + |Au,„p, 

then by ^ 

(97) ^ <C(uo)2/'^+^ 
Integrating this differential inequality, we find 

0<y(t)<(2/(0)-W)_c(uo)t)-^/('''+'^ 
for < i < To where 

< To < ^ 



C(wo)(l + |Auo|2)4(p+i)' 

This together with (|M)) implies that as in (jTO)) satisfies > Tq for each m 
and 

(98) Mm G a bounded set of ^^(0, Tq; W'l) n i°°(0, Tq; M^), 

independent of m. 



Now, by (|7^ and ((M)) we have the following estimate for |^ff-|i2(o 



(99) <C(^io)(|A2w„||2(o,To;L2) + |VUm • V Aw,„ | ^2 (q, Tq ; L^) 

+ \A{biUra + b^ul^ + &3"m)li2(o^To;L2) 
+ |Vli„,. • V(6lM,„ + b2U^n + &3Mm)|i2(o,To;L2))- 

Note that 

|VUm-VAUm|i2 < IWrnl/i-i <C| W,„ | l | Mm | 

^C'(|Wml/i"i + l^mlffi); 

which together with ([83l) implies 

|VUm • VAu™|2 2(o_y^.^2) < CiC{uo)To + |A2u„|2,^Q^^^.^,^). 
Similarly, we have 

|A(6ilim+62W^, + &3Wm)li2(o,To;L2) < C(C(uo)To+ | A^ W„ | ^2 (q, Tq; L^) ) ' 

and 

|VUm • V(6ili,„ + 62Mm + Mm) |i2(o, To; L^) < C(C(uo)To + | A^Wm |i2 (q, To; L^) ) ■ 



26 



LIU, SENGUL, AND WANG 



Plugging the above four inequahties in (|99|l . we obtain 



(101) 



which together with (|M|) shows that 

(17/ 

(100) ^ ^ ^ bounded set of L^O, Tq; i^). 

Step 2. By (1^ and (jlOOp we can extract a subsequence Um' of Um which 
satisfies 

Um' u weakly in L^(0, Tq; Wi), 
Um' u weak-star in L°°(0, To; W), 

^ ^ ^ weakly in L'{Q,To\L''). 
At dt ^ \ ^ V, J 

Thanks to the compactness of the embedding of Wi in H'^ f) W, the inclusion 

{f e L\0,To;W,) I % e L\0,To;L^)} C L\0,To;W n H^) 
at 

is compact; see e.g. [3] Lemma 1.6. Therefore without loss of generality, we may 
assume 

(102) Um' ^ u strongly in L^{0, To; W n H^). 
Note also the following embedding is continuous 

{ / e L^O, To; W^i), ^ e L\0, Tq; L^) } ^ C([0, Tq]; W), 
see e.g. [3]. Thus, upon passing to a further subsequence, we have by poip i and 

mils 

(103) Um' u weakly in C([0, To];W). 

In particular, Mm'(O) converges weakly to u(0) in VK, and so u(0) = uq because 
Wm'(O) converges to uq strongly in W. We still need to show that the function u 
satisfies ((M|) 9. 

We consider cj) e C;f (0,To) and TV > 1. For any m' > TV, satsfies ^2 with 
w = cn where e^r is as in (j90p . We multiply this equation by </)(<) and integrate by 
parts to obtain 

(104) 

Um'EN'P' dt 

; 

= / {H{Um')f^m'^eN + H'{Um')^J.m''^Um' ■yeN)4ldxdt. 

JO Jn 

The convergence properties of the sequence Um' allow us to pass to the limit in this 
equation. The passage to the limit on the LHS is easy to see by using (|10ip i . and 
we have 

(105) / / Um'eN<P'dxdt > / / mcatc/)' dx dt. 

Jo Jn Jo Jn 
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For the RHS, we have 
(106) 

To 



"'SI 

To 

a 



Jo Jn 



10 Jn 

To p 

/ H{Ujn'){blU,n' +&2Wm' + 63?/^,, ) Abat^ dx dt 

Jn 
Jn 

For brevity, we wiU only show the convergence of the first term and the convergence 
of the rest terms fohows m the same fashion. 

pTo p pTo p 

I / / H{um')Au„i' AeM<pdx dt — / / H{u)AuAeM4'dxdt\ 
Jo Jn Jo Jn 

(107) ^ I / " / {H{u„i') ' H{u))AuAeN(j)dxdt\ 

Jo Jn 

To p 

/ H{u){Au„i' - Au)AeN(tidxdt\. 
Jn 

Using (jlOip - (|103p and mean value theorem, the first quantity on the RHS of (|107p 
can be estimated as 



(108) 

I / / {H{u„i') - H{u))AuAeN(f>dxdt\ 
Jo Jn 



< I / / H'{wm'){um' — u)AuAeN<pdxdt\ 
Jo Jn 

<C [ [ Qwm'lP + l)\u,n' - u\\Au\\AeN\\(t)\dxdt 
Jo Jn 

< C(|u;™'|^„(Q,J,^.^y) + l)|u|io=(o,To;VK)l0|L°° J \ujn' - u\\AeN\ dx dt 

< <^(kmnL-(0,To;H^) + l) I '^U- (0,To 1 0| / l^m' " uIl^IcnIh^ dt 

J 

< C(|Wm'lLoo(o,To;W) + 1) I "I L~ (0,To ; W) 1 (^l^m' - "lis + (5|eAr | ^2 ) dt 

< C'<5(|w„i'|^^(.g y^.^y^ + l)l"lL==(0,To;W)l'/'l-L°°|ew|i~(o,To;W) 

c 

+ (l^™'li=o(o,To;W) + l)l'"lL~(0,To;H')l</'|L~|Wm' - ^1^2 (o,To;W) ■ 

In light of (|102p . the above quantity can be made as small as possible by choosing 
S > sufficiently small and to' sufficiently large. The second quantity on the RHS 
of (|107p can be estimated in the same way. 
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Now, we obtain after passing to the limit the fohowing equation for u: 

— I I ueN(t>' '^x dt = / / {H{u)ijAeN + H' {u)ii\7u ■ VeN^ipdx dt. 
Jo Jn Jo Jn 

The hmit equation obtained above is fulfihed for any N and any (f) G C^{0, Tq), so 

that the density of span{eN \ N € N} in W allows us to conclude that u satisfies 

(1111)2. 

Step 3. In the following, we will sketch the proof for the uniqueness. Let ui 
and U2 be any two strong solutions of (|74p defined on the interval [0, Tq]. There 
exists C(To) > such that for z = 1, 2 

(109) \uS)\h- <CiTo), VtG[0, To]. 

Let u ~ ui — U2- Multiplying (fM)) by w G i/^, integrating over £7, we get: 

(110) v) ={H{u)\/{aAu - {hu + b^u^ + 63"')), Vw), 
dt 

from which we see that u satisfies 
du 

(— , v) ^{aH{ui)WAu, Vv) + {a{H{ui) - H{u2))VAu2, Vv) 

•^m) - {H(ui){biUi + b2u\ + bsul - biU2 - b2ul - 63U2), Vw) 

- {{H{ui) - H{u2)){biU2 + b2ul + 63^^), Vf). 

From the regularity we obtained for solutions of (|74|) . we can take v in the above 
equation to be —Au and use {H) to get: 



1 dlVw 



|2 



~aBi\VAu\^ < -~{a{H{ui) - H{u2))VAu2, VAu) 



2 dt 

(-^^2) + {H{ui){biui + 62-1*1 + &3M? - biU2 - b2ul - b^ul), VAu) 

+ {{H{ui) - H{u2)){biU2 + b2ul + b?.ul), VAu). 

Here, the term '^^^"^ is understood in the distribution sense. More specifically, 
since ^ £ L^{0, To;L'^) and u G L^{0, To;Wi), then by Theorem 2.3 in [13], we 
know that ^ G L'\0,To] W). 

Denote the terms on the RHS of (|112[) by 13, I4, I^. We have the following 
estimates for them. 

Applying mean value theorem to H and using (|109p . we have 

/3 = - {a{H{ui) - H{u2))VAu2, VAu) 
<C\H'{w)\l'^ |u|loc IVAU2I |VA{t| 
<(by Agmon's inequality) 

<C(1 + \w\l^)\u\]j, \u\l, IVAU2I |VAm| 

<C(1 + \w\l^)\ii\l, \u\l, \VAu2\ \VAu\ 

<C{1 + \wf^^,)l\\7Au2\^ \u\'hi + ^iVA^ip 



<c(i + c(Tor)H\u2\i. + mil. + ^\^Au\ 



where w = 6{t)ui + (1 — 0{t))u2 for some 0{t). 
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By UMl), wc have 

li ^{H{ui){hiui + b2ul + - hiU2 - - ^3^2)- VAu) 
<C(1 + W\ltt){l + \u,\l^ + |^.2|ioo)|S||VA^| 

<C{1 + C{ToY+^f\u\^ + ^|VAu|2 



Similarly, 



<C(1 + C(To)P+3)2|V?ip + ^iVAup. 

/g ={{H{ui) - H{u2))ibiU2 + b2uj + bsul), VAu) 
<C(l + C(To)P+3)|u||VAu| 

<C(1 + C(Tof+3)2|Vu|2 + ^|VAu|2. 
Plugging the above estimates in (|112p . we have 

(113) ^^<Ci\u2\l. + l)\Vuf, 

which together with u € L'^{0, Tq; Wi) and |Vu(0)|2 = implies \Vu(t)\'^ = for 
all t € [0, To], and the uniqueness is thus proven. 

□ 

Completion of the proof of Theorem \ 6. 1[ For 1 < p < 3, one can find 2 < 
92 < 3 such that the following inequalities are satisfied: 

(114) -"-"^ ^2 ^ 

3(p + 3)<-^, b + 3)(|-l)<2. 
3-92 V 2 92 / 

Taking inner product on both sides of ((74)) with A^u, we get: 

i^lAup = - {aH{u)A^u, A^u) - {aH'{u)Vu ■ VAw, A^u) 

(115) + {H{u)A{biu + b2u^ + bau^), A^u) 

+ {H'{u)Vu ■ V(6iu + b2U^ + bsu^), A^u). 

Here i^|Aup on the LHS is understood in the scalar distribution sense on (0, T); 
again see Theorem 2.3 in |13j . Then by our assumption (H), we have 

i-^|Au|2 + a^ilA^ul^ < -{aH'{u)Vu ■ V Au, A^u) 

(116) + {H{u)A{biu + b2U^ + 63"^), A^u) 

+ {H'{u)Vu ■ V{biu + b2u'^ + 63"^), A^u), 
:= Ji + J2 + ^3- 

Let pi = 6p||±||. Then for ^ > sufhciently small we have pi < ^ by pii)) . 
Hence l^^^'"?! '-^ L^^ and we have 



(117) < C\Vu\Ln < C\Vu\'^~^\yu\l, " < C(uo)|u||, " 
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Using ()117|) and ()114p we can obtain: 



<C{un)il + \u\j,,). 
To estimate Ji, let 77 be defined as 



3 _ 3 



3 + 13 6 - ?/ 2 
and note that 

(120) |VAm|^«!', < ivAT/r/^ivAT/i^;'"/" < 

We estimate Ji using (jllSp . (|119p and (|120p as follows: 
Ji = - (aff'(w)Vu • VAm, A^u) 
<C|i7'(u)Vw|i3+3 |VAu|i6-„ I A^ul 

<CK) (1 + \u\l.) (,-(12(6-.)-.)/.) + ,|^|2^^)_ 

By pTi)) . wc know M^^'^^ ^3(p+3)^ ^j^g^^ 

(122) |mL3(p+3) < C\Vu\l.2 < C\Vu\^^^\Vu\i;^ < C{uo)i\uf^^_P+^^ - 
To estimate J2 we first estimate the following two integrals 

/ (1 + \u\P+')\Au\ \A\\ dx <C{1 + |«|^tf,^3,)|A«k« lA^^I 

<by (inSI) 

(123) <C{uo){l + \u\j,,)\u\H^\u\H^ 

<Ciuo){l + \u\l.)\u\]/f\u\]/! 
<C{uo){l + \ufH.){e-' + e\u\l.) 

f {l + \u\P+^)\Vuf\A^u\dx<C f (l + |u|P+3)|VunA2u| 
Jn Jn 

<C{l + \u\lZ^,,)\Vu\UA^u\ 

(^2^) <Ciuo){l + \u\l,)\u\l,Mu\H^ 

<Ciuo){l + \u\l.)\uf^^\u\'^{' 
<C{uo){l + \u\l,){e-''+e\u\l,). 



Using p23)) and p24)) we have 

J2 = {H{u)A{biu + b2U^ + bsu^), A^u) 

(125) <C / (1 + lul^^+i) [(1 + |up)|Au| + (1 + |ti|)|Vu|2] lA^ul dx 
Jo 

<c(«o)(i + Hl^O(e"" + e|«l?/4). 
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For J3, we have 

J3 ^{H'{u)Vu ■ V{biu + b2U^ + hu^), A^u) 

< J (l + |u|P)(l + |u|2)|Vu|2|A2u|dx 

<(by ) 

<Ciuo){l + \u\l.){e-''+e\u\l,). 



(126) 



From the estimates for Ji, J2 and J3 given in (|12ip . (|125p and (jl26p respcetively, 
we have by ^ and pTC)) : 

(127) - C(wo)e^(i)) < C{uo)e-^ A{t). 

Here = maxjll, (12(6 — f]) — rf)/ri} with determined by ()119p . Also note that 
N ^ 00 as p approaches the critical exponent 3. 

The crucial step towards the global existence and uniqueness is a uniform 
bound for the solution. This can be achieved by manipulating (jl27p when the initial 
data is small as we now show. To our knowledge, a similar method first appeared 
in [12]. 

First, for any t>Q and 1 > e > to be specified later, we have by ()83p and ([Si)) : 

/ ^ A{t) dr = /" ^ {\u\]j, + 1) dr < {C\u\hi \u\h^ + 1) dr 
Jt Jt Jt 

(128) </ iC\u\jj,+C\u\jjs+l)dT 

<Ce{l + \u0\l2f + Cluoll^il + \u0\l2) 



+ ci{i + \uo\i,2r' + 



e. 



From now on we will assume that |wo|//2 < 1. Then we have by (jl28p 
(129) A{T)dT<C{i+\u„\l.), 



where C is independent of uq. 
Let 



Ci=C{uo)e, C2^C{uo)e- , 
(130) 

e~=e'^, M^C{i+\uo\l2). 

It is easy to see that there exists e > sufficiently small such that for any initial 
data uq satisfying |mo|^2 ^ we have 

(131) aBi > Ci(^(0) +C2A/ + 4C). 

Then by the local well-posedness, we know that there exists T* > such that 

(132) aBi > CiAit), for t < T*. 

We claim that T* > e. Otherwise, by (fT27)) . (fT30)) and (fT32)) . we have 

(133) ^<C2^(t) V<G[0, n, 

at 
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then by ^29]) 

(^34) A{T*) ~ ^(0) < C2 1^ Air) dr <C2 |' ^(r) dr 

which leads to the following contradiction to the definition of T* : 

(135) aBi > CyA{T*). 

We claim now that T* = oo. Otherwise, by 3 i* £ [T* - |, T*], such that 

(136) A(t*) < AC. 
We also know 

(137) A{T*) - A{t*) < C2M. 
Thus 

(138) A{T*) <4C + C2M. 

Again, we are led to the contradiction (|135p . 
Since T* = cxd, then 

(139) aBi > CiA(t) ^ Ci{\u{t)\jj2 + 1) Vt>0, 

which implies the uniform bound of the solution. 

Finally, Theorem 16.11 follows from Proposition 16.31 and (|139[) . □ 

6.4. Proof of Theorem 16.21 We first give a lemma on the existence of solutions 
to the following Cauchy problem; 

(140) dt ^ ■/ V 7, 

w(0) = uo, 

where A is the infinitesimal generator of an analytic semigroup in X with domain 
D{A). 

Lemma 6.3. Let uja = sup{ReA | A G (j{A) } < 0, / G Cfc([0, 00); Da{0)), and 
uq G Da{0 + 1), where (j{A) is the spectral set of A, DA{d) and Da{0 + 1) are as 
defined in (jl4p with some Q < 9 < 1. Then there is a unique solution of (|140[) which 
belongs to Cfc([0, cxd); Da{Q + 1)) H C^([0, 00); Da{0)), and there exists a constant 
C independent of f and uq, such that 

^^^^^ II"IIc6([0,oo);_Da(9+1)) + \W\\cMO,<x,);Da{0)) 

< C'(ll/IICi,([0,oo);DA(e)) + \\uo\\da{0+1))- 

This lemma is a direct consequence of sections 4.3 and 4.4 of Lunardi [15]. 

We first show that Theorem 16.21 is true when T > Tc- In this case, from ((37)) 
we see that ujl^ = sup{ A | A G cr{LT) } < 0. Now for any given v G B{0, 1) C 
C(,([0, cxd); Dl^{9 + 1)), we consider the following linear equation: 

(142, |^L.. + G(„,r), 

u(0) = mq. 
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With our choice of 6', the space Dlj.{(^) forms an algebra according to Lemma 
14.11 Then it is easy to see that G{v,T) e Cb{[0, oo); 1)^^(0)). For example, the 
term v{t)'^ A'^v{t) in G(u,T) can be estimated as 

\\v{t)'AMt)\\D,^ie) < ll«Wlb,,wlk(i)lb,,wl|A't'(i)lb,,w 
< \Ht)\\l,^ieMt)\\D,,ie+i) 

So u^A^w £ Cb([0, oo); 1)^^(6')). Applying similar estimates to other terms in 
G{v, T) we obtain the following: 

l|G'(w,r)||c,([o,oo);r>i^(e)) <Cilkllci,([o,oo);_Di^(e+i)) + o(||w|lc,([o,oo);Di,^(e+i)))' 

where Ci is independent of v. 

Now, by Lemma lOl (|142p has a unique solution m in C(,([0, oo); Dl^{9 + 1)) n 
C^([0, oo); Diy(6')), which satisfies 

ll^i|Icb([o,oo);Dtj,(o+i)) <C' { ||G(u,T)||c^([o,oo);Di,^(e)) + ll"o||r>i^(e+i) ) 
(143) ^ ^ 

<C2 (lkllc,([0,oo);Di^(e+l)) + ho||Di^(9+l) j • 

Let i? = min{^^, 1}, Bi be the ball centered at zero in Cb([0, oo); Dl^{6 + 1)) 
with radius R , and B2 be the ball of radius centered at zero in Dl^{9 + 1). 
Define a mapping F as follows 

F:BixB2^Bi, T{v,uq)=u, 

where u is the solution of (|142p with given v and ug. By (|143p and our choice of 
i?. F is well defined. 

Now we will prove that F is a contraction in the first variable. For any vi ,V2 S i?i , 
let T(vi,uo) = Ui, i 1, 2. Let u = ui — U2 and v ~ vi — V2. Then u satisfies the 
following equation: 

(111 

(^44) -^=LTU + Giv,,T)-Giv2,n 

u{0) = 0. 
Again by Lemma l6.31 we have 

ll"llc6([o,oo);_Dt^(e+i)) <C\\G{vi,T) - G{v2,T)\\c,([o,oc>}-DL^{e)) 

<C2\\v\\c,{[o, 00); DLj.{e+i)){\\vi\\cMo., 00); _Di.^(e+i)) 

+ l|l'2||c,([0,oo);_Di^(e+l))) 
<2i?C2||w|jc5([0,oo);_Dt^(e+l)) 

1.. I, 

Namely, 

||F(wi,Mo) - r(w2,Uo)||c,([0,oo);_Dz,^(e+l)) < ^11^1 " "2 || Ct,([0. oo) ; Dz,^ (6+1)) ■ 

From above, we see that given any uq G B2 there is a unique fixed point u G _Bi 
such that F(w, tin) = w. So for any initial datum uo £ B2 C Dlt{0 + 1)7 the 
equation (j80| admits a unique solution u E C(,([0, 00); Dlj.{& + 1))- It is easy to 
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see that u is also in C^([0, oo); Dlt.{S))- The theorem is proved in this case with 
r = R^. 

For the case when T <Tc, wc define 

Lt ^ Lt - iPiiT) + 6) id, 

G{;T) = G{;T) + {Pi{T) + S)id, 

where id is the identity map, (3i{T) is the largest eigenvalue of Lt, and 6 is some 
positive number to be chosen below. 

By p7)) . we can choose e and S sufficiently small such that 

(145) \(3iiT) + S\ < R, y T e [T, - e,T,], 

where R = minj^j^, 1} as before. 

Now consider (|142p with Lt and G replaced by Lt and G, respectively. Note 
that w^;^ = sup{ A | A £ <j{Lt) } < 0. Following the same argument as for the case 
T > Tc with suitable modification and making use of (|145p . one can show that for 
any uq G B{0,R'^) C Dlt{0 + 1), there is a unique u 6 Gbi[0, oo); Dlt{0 + 1)) n 
GliiO, oo); Dlj.{0)) such that 

^ =Ltu + G{u,T) 
at 

^Ltu - (/3i (T) + 5)u + G{u, T) + (/3i (T) + 5)u 
=Ltu + G{u,T). 

The proof is now complete. 

References 

[1] J. W. Barrett and J. F. Blowey. Finite clement approximation of the Cahn-Hilliard equation 
with concentration dependent mobility. Math. Comput., 68:487-517, 1999. 

[2] Peter W. Bates, Kening Lu, and Chongchun Zeng. Existence and persistence of invariant 
manifolds for semiflows in Banach space. Mem. Amer. Math. Soc, 135(645);viii+129, 1998. 

[3] F. Boyer. Mathematical study of multi-phase flow under shear through order parameter for- 
mulation. Asymptot. Anal, 20(2):175-212, 1999. 

[4] J. Cahn and J. E. Hilliard. Free energy of a nonuniform system. I. Interfacial free energy. J. 
Chem. Phys., 28:258-267, 1958. 

[5] G. Da Prato and P. Grisvard. Equations d'evolution abstraites non lincaircs de type 
paraboliquc. Ann. Mat. Pura AppL, 120:329-396, 1979. 

[6] G. Da Prato and A. Lunardi. Stability, instability and center manifold theorem for fully 
nonlinear autonomous parabolic equations in banach space. Arch. Rational Meek. Anal., 
101:115-141, 1988. 

[7] H. Dang. Stability and boundary layer properties of Cahn-Hilliard equations. Ph.D. Thesis. 

University of Utah, 1995. 
[8] C. M. Elliott and H. Garcke. On the Cahn-Hilliard equation with degenerate mobility. SIAM 

J. Math. Anal, 27:404-423, 1996. 
[9] C.M. Elliott. The Cahn-Hilliard model for the kinetics of phase separation, in Mathematical 

Models for Phase Change Problems, Internat. Ser. Numer. Math. 88, Birkhduser-Verlag, 

Basel, Switzerland, pages 35-73, 1989. 
[10] D. Henry. Geometric theory of semilinear parabolic equations, volume 840 of Lecture Notes 

in Mathematics. Springer- Verlag, Berlin, 1981. 
[11] C. Hsia. Bifurcation of binary systems with the onsager mobility. J. Math. Phys., 51:063305, 

2010. 

[12] F. H. Lin and C. Liu. Nonparabolic dissipative system modeling the flow of liquid crystals. 
Commun. Pure Appl. Math., XLVnL501-537, 1995. 



CAHN-HILLIARD EQUATION WITH ONSAGER MOBILITY 



35 



[13] J.L. Lions and E. Magenes. Nonhomogenous Boundary Value Problems and Applications. 
Springer, Berlin, 1972. 

[14] C. Liu, Y. Qi, and J Yin. Regularity of solutions of the Cahn-Hilliard equation with non- 
constant mobility. Acta Mathematica Sinica, English Series, 22(4):1139— 1150, 2006. 

[15] A. Lunardi. Analytic semigroups and optimal regularity in parabolic problems. Birkhauser, 
1995. 

[16] T. Ma and S. Wang. Phase transition dynamics in nonlinear sciences. To appear. 

[17] T. Ma and S. Wang. Bifurcation theory and applications, volume 53 of World Scientific Series 
on Nonlinear Science. Series A: Monographs and Treatises. World Scientific Publishing Co. 
Pte. Ltd., Hackensack, NJ, 2005. 

[18] T. Ma and S. Wang. Cahn-Hilliard equations and phase transition dynamics for binary sys- 
tems. Disc. Cont. Dyn. Sys. B, 11:741-784, 2009. 

[19] T. Ma and S. Wang. Phase separation of binary systems. Physica A, 388:4811-4817, 2009. 

[20] A. Mielke. Locally invariant manifolds for quasilincar parabolic equations. Rocky Mountain 
J. Math., 21(2):707-714, 1991. 

[21] A. Miranville and S. Zelik. Exponential attractors for the Cahn-Hilliard equation with dy- 
namic boundary conditions. Math. Methods Appl. Sci., 28:709-735, 2005. 

[22] A. Novick-Cohen and L. A. Scgcl. Nonlinear aspects of the Cahn-Hilliard equation. Physica 
D, 10:277-298, 1984. 

[23] G. Schimperna. Global attractors for Cahn-Hilliard equations with nonconstant mobility. 

Nonlinearity, 20:2365-2387, 2007. 
[24] G. Simonett. Center manifolds for quasilincar reaction-diffusion systems. Differential Integral 

Equations, 8:753-796, 1995. 
[25] R. Temam. Infinite- dimensional dynamical systems in mechanics and physics, volume 68 of 

Applied Mathematical Sciences, second edition. Springer- Verlag, New York, 1997. 
[26] H. Triebel. Interpolation Theory, Function Spaces, Differential Operators. North Holland, 

Amsterdam, 1978. 

(HL) Department of Mathematics, Indiana University, Bloomington, IN 47405 
E-mail address: liu409indiana.edu 

(TS) Department of Mathematics, Indiana University, Bloomington, IN 47405 
E-mail address: msengul9indiana.edu 

(SW) Department of Mathematics, Indiana University, Bloomington, IN 47405 
E-mail address: showaiig9indiana.edu, http://www.indiaiia.edu/_fluid 



